Equation of motion in a scalar model of 

gravity 

Q^ ' Shmuel Kaniel and TYakov Itin 

^_» ' Institute of Mathematics 

o : 

^^ , Hebrew University of Jerusalem 

1/-^ . Givat Ram, Jerusalem 91904, Israel 

^ ; February 7, 2008 

>■ 
m ■ 

00 . 

O. 

O ' Abstract 

y—i ] 

^^ ' A scalar model of gravity is considered. We propose Lorentz invariant field equation 

— , . □ / = kr]abf,af,b- The aim of this model is to get, approximately, Newton's law of 

^^1 gravity. It is shown that / = — ^ln(l — k^) is the unique spherical symmetric static 

' ■ solution of the field equation. / is taken to be the field of a particle at the origin, 

bj). having the mass m. The field of a particle moving with a constant velocity is taken 

^1 to be the appropriate Lorentz transformation of /. The field F of N particles moving 

IT^ ' on trajectories ipj{t) is taken to be, to first order, the superposition of the fields of the 

'j_j ■ particles, where the instantaneous Lorentz transformation of the fields pertaining to 

. . . . the j-th particle is ipj{t). When this field is inserted to the field equation the outcome 

is singular at {ipj{t),t). The singular terms of the l.h.s. and of the r.h.s. are both 

0{R~^). The only way to reduce the singularity in the field equation is by postulating 

Newton's law of force. 

It is hoped that this model will be generalized to system of equations that are covariant 
under general diffeomorphism. 



1 Introduction 

Classical field theory postulates the world to be a 4D-differential manifold M. The evolution 
in M is represented by field equations. The most fundamental are: the scalar wave equation. 
Maxwell and Einstein field equations. To these one may add Newton equation for the 
gravitational potential. It is not an evolution equation, but a field equation for a section 
t = const. For the field equations above particles are postulated to be the singularities of 
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the field. Newton and Maxwell equations are linear. Therefore, there is no way to formulate, 
only by the field equations, the "force" that a particle exerts on other particles. Thus 
additional laws were included: Newton's law of gravitation. Coulomb and Lorentz laws in 
electrodynamics. 

On the other hand, Einstein field equation for gravitation is non linear. This allows, in 
theory, to derive equations of motion from the field equations. Einstein succeeded to carry 
out this program. The original version of general relativity (GR) employed an additional 
postulate: The Geodesic Postulate. Thus, he obtained a covariant generalization of the 
classical equation of motion. The Geodesic Postulate is an intrinsic postulate involving only 
the geometry of M. By comparison, Newton, Coulomb and Lorentz laws are extraneous laws 
of force. 

Later, Einstein- Infeld-Hoffman [Q] and Fock [^ derived, directly, the equations of motion 
from the field equations. The derivation is rather formal (i.e. the first terms of a series). 
Moreover, it is not known how the motion of the particles can be embedded in a field 
satisfying Einstein field equations. Recently Damour [Q] computed the equations of motion 
to great accuracy. Still, the derivation is formal. 

In this paper we study a Lorentz-invariant scalar model of gravity. This is a non-linear 
(quadratic) Lorentz-invariant generalization of the Newtonian scalar theory of gravity. We 
derive by this model the equations of particle motion and Newton's law of force from the 
field equation. 

Recently, Watt and Misner ||K| considered a scalar model of gravity (It is interesting to 
note that the metric that they obtained is, also, the metric obtained by Kaniel and Itin |P 
i.e. Yilmaz-Rosen metric.) The motive of [^] was to facilitate numerical computations of 
gravitational waves. We hope that the model established in this paper will be of case also for 
analytical reasoning. It may serve as guideline for studying the equations of motion, where 
the field equation is a diffeomorphic covariant system. We propose to apply the method 
exhibited in this paper to general quasi-linear field equations. The computations above lead 
to a novel algorithm for the derivation of equations of motion from it. 

1 Compute a static, spherically symmetric solution of the field equation. It will be 
singular at the origin. This will be taken to be the field generated by a single particle. 



2 Move the solution on a trajectory ip{t) and apply the instantaneous Lorentz transfor- 
mation based on '0(t). 

3 Take the field generated by n particles to be the superposition of the fields generated 
by the single particles. 

4 Compute the leading part of the equation. Hopefully, only terms that involves ifj will 
be dominant. 

5 Compute the "force" between the particles by the quadratic part of the equation. 

This algorithm should result, in most cases, in the verification of Newton gravitation law. 
The method presented in this paper may enable the embedding the trajectories satisfying 
the equations of motion in a field satisfying the field equations. 



2 Non-interacting particles 

We start with the hnear d'Alembertian field equation. This equation is a Lorentz-invariant 
extension of Laplace equation pertaining to the Newtonian gravity. 

The Newtonian theory of gravity is formulated as a field theory via the scalar potential /, 
which is subject to Laplace equation 

A/ = 0. (2.1) 

The unique spherical symmetric and asymptotically zero solution of this equation is 

/ = i. (2.2) 



where R = va^^^oP) ''^ ^^^ ^o are arbitrary constants of the integration. This solution 

is singular, so it represents a field of a pointwise particle located in the point f = Tq. The 

parameter m, taken to be positive, represents the mass of the particle. 

Laplace equation is not Lorentz-invariant. Thus it should be generalized to the wave field 

equation 

n/ = 0. (2.3) 

The time independent spherical symmetric asymptotically zero solution of this equation is 
again 

f-^. (2,4) 



where 



Jl= (^r -ro) -v(a <v,{r - tq) > +/3tj , (2.5) 



provided that ( |2.5| is a Lorentz transformation (a and (3 are functions of v). 

We reaffirm this fact in order to establish the notation. 

This solution describes free motion of a particle. Indeed, the singularity point satisfies 

f — fo = via <v,{f — fo) > +pt] . (2.6) 

So 

<v,{r- To) >= v^ ia <v,(f- Fq) > +/3t j . (2.7) 

Or 



1 — av'^ 



< ^;, (r - ro) >= ^ j7^. (2.^ 



Substituting this relation in (|2.6|) we obtain the equation motion of the singularity as 



/3 ^. 
1 — av^ 



f-fo = -r^^vt. (2.9) 



In order to describe by ( [^.9| ) a constant velocity motion we have to take 

/3=l-at;2. (2.10) 

The functions a, (3 are still arbitrary. However the field equation (|2.3| has to be satisfied. 



Calculate the derivatives 

Rt = -v(3, Rtt = 0, (2.11) 

Rx = ei- aviv, R^x = 0. (2-12) 

Thus 

ft = -^<Rt,R>=^<v,R> (2.13) 



and consequently 



fjlQ TfiB 

ftt = -^ <v,Rt> -S^ <v,R>< Rt,R> 

'^^'^3<v,R>^-v^R'). (2.14) 



i?5 

As for the spatial part of the d'Alembertian 

fx = -^<Rx,R>, (2.15) 

fxx = ~~m ^ -^2;, Rx > +3-^5 ^ -^2;, R > ■ (2-16) 

Thus the Laplacian 

Af = 3^(^<R,,R>^ + < Ry,R >2 + < R,,R >^^ ^ ^(^Rl + Rl + R^}j. (2.17) 

Using the relation 

<Rx,'R>= Ri-avi<v,R> (2.18) 

we obtain 

<Rx,R>'^ + <Ry,R>^ + <R„R>^=R^-2a<v,R>^+a\^ <v,R>^ (2.19) 
and 

Rl + Rl + Rl = 3- 2av^ + a^v^. (2.20) 



Hence the Laplacian takes the form 

R^ it"* 

= ^{a^v'^ -2a)(?,<v,R>'^ -v'^R^\ (2.21) 

Since D / = 0, equations ( p.l4| ) and ( |2.21| ) imply that 

/?2 = q;2^2 _ 2^_ (2.22) 

By ( p.l4| ) the Lorentz parameters a and (3 are 



/3 = ^=i=^, « = 4 (l - ^=^) (2-23) 

The result is rather obvious. Indeed, beginning with the static 1-particle solution ( p.2| ) of 
the wave equation ( p.3|) and making a Lorentz transformation of coordinates (with opposite 
velocity) we obtain a solution which describes the inertial motion of the particle. 
Consider now the function 

/ = 9, (2.24) 

where 

i? = (r - To) - a-0 < -0, (r - Fq) > -(3^1, (2.25) 



with il) = il)[t) and a and /? are functions of 



12. 



/?=^i^, « = ^(1- ^-^^ ) (2.26) 



i-hAp i^r a/1- 



If ■?/' = f t then (p.5|) coincides with 

Calculate the d'Alembertian of the function /. 

Rt = —2a' ( < ■?/', -0 > i/j < ■0, (r — fo) > j — a ( -0 < t/^, (r — Fq) > 

-a('^<^,(r-ro) >) - 2/3'(' < ^,^ > ^) - /3^ (2.27) 

and 

Rtt = -4a"(<^,^>2<^,(r-ro) >) -2a'(|^|2<^,(f-fo)>) 

-2aM < t/;,?/) >< i/j,(f -To) > j - 2a (< ip,ilj >< i/j,(f -Tq) > j 

-4/?"('<^,^>2^) -2P'(\^\'^^) -2/5'('<^,^>^) -2/3'('<^,^>^) 
-2/?' f < ^, ^ > ^) - /3^. (2.28) 



Thus 

fjl JTi ( /J— .11- J— _:_ 

ft = -^<Rt,R>=-^{ -2a i^<'i/),ilj><ilj,R><'ilj,{r-ro)> 
a( <iIj,R >< ilj,(f -To) >j - a( <iIj,R >< i{j, (r - rg) > ) 
-2P'(<i^,i><i^,R>) -l3(<^,R>)] (2.29) 

In the first approximation (for tlie motion of tlie particle, wliicli is slow with respect to the 
speed of the light) we can take 

Rt = -P'^, Rtt = -13^. (2.30) 

Consequently, 

ft = ^P<i^,R> (2.31) 

In general, the derivatives of a and [3 contribute terms that are quadratic in if) and its 
derivatives and can be neglected. Thus the second derivative to the same accuracy is 

fflR / .11. _:_ \ TTlB -^ 

ftt = -^[<^,R> + <i',Rt>)-S^<iP,R><Rt,R> (2.32) 

Substitute 

Rt = -0^ (2.33) 

to get 

ftt = -^[<^,R>-P <^,^>) + 3^<ij,R>^ 

= ^<tR> + ^{3<i^,R>'-R'\W)- (2.34) 



As for the spatial derivatives 

Rx = ei- aipii/j, Rxx = (2.35) 

777, 

fx = -j^<Rx,R> (2.36) 

/rn 777 

fxx = ~"m ^ ^^' Rx > + 3— < Rx, R > (2.37) 

A f = 3^(^<Rx,R >^ + <Ry,R >^ + <Rz,R >^ ) - ^[rI + Rl + Rl) (2.38) 



Since 

<R^,R>=Ri-atpi<t,R> (2.39) 

we get 

<R^,R>'^ + <Ry,R>^ + <R„R>'^=R'^-2a <^,R>^ +a^\^\^ <^,R>^ (2.40) 
and 

Rl + Rl + Rl = 3 - 2a\^\^ + a^\ip\' (2.41) 

Thus the Laplacian is 

A/ = ^,{aW - 2«) (3 < ?, /? >' -m'R') (2.42) 

And, for the d'Alembertian 

□/ = ^ < ^, i? > +^iP' + 2a- a^m (3 < ^,i? >' -R'\W) (2.43) 

Using the expressions for the functions a, P we obtain 

Df = ^<t,R> (2.44) 

Thus in the case of a hnear field equation D/ = the motion of a one singularity point is 
inertial i.e. ^ = 0. 

The reasoning above can be extended to arbitrary number of singularities. Indeed, the 
equation ( |2.3| ) is linear so it has a solution 

1=1 
where 

(^R = (r - (*)r) - »a(^)^ < (^)^, (f - (^V) > -(^)/?»^, (2.46) 

with 

«/3 = ^i=, (^^a=^(l- ' . ) (2.47) 

The calculation above yields, approximately 

n „ 

n/ = E^<^^'^^>- (2-48) 

2=1 * 

Also here, for the field equation is D/ = / is a solution provides that -0^ = 0. 
Thus, the linear equation describes the free inertial motion of an arbitrary system of singu- 
larities. 



3 Non-linear equations 

In order to describe the interaction between singularities we need to balance the value ( |2.48| ). 
Let us consider a new non- linear field equation 

□/ = kv^'laU (3.1) 

where /c is a dimensionless constant. It is easy to see that this equation is unique Lorentz- 
invariant equation which is linear in the second derivatives and quadratic in the first deriva- 
tives. 

The first step is to find a static spherical-symmetric solution of ( p.l|) . Write / = f{s) with 
s = x"^ + y'^ + z'^) . So 

Af = 2f + Af"s (3.2) 

and 

V'^'laU = -4(/')'s. (3.3) 

Thus the equation takes the form 

2fs + 3f = 2kifys (3.4) 

Inserting f = Z we obtain 

2Z's + 3Z -2kZ^s = (3.5) 

Take the new variable 

Z = s"^/2y (3.6) 

to obtain 

Y' = ks-^'^Y\ (3.7) 

Thus 

Y = ^^ ===> f = (3 H) 

C^s + 2k ■' s{Cy^ + 2k) ^ • ' 

And the unique exact solution of (|3.1|) is 

/--j'"(i-^v)- l^'S) 

Note that for the limit A; — > we obtain / — > ^. The solution ( p.9| ) is singular at the point 
r = 0. The singularity can also be located in an arbitrary point r = tq. 
The next steep is to consider a moving singularity. 
We seek for a solution of the form 



for a moving singularity: 

R = (f — fo) — ai/j < -0, (r — Fq) > — /3'0, (3-11) 

where a and (3 are functions of j-j/'p. From the calculations above 

Rt = -P'i, Rtt = -0i. (3.12) 

Thus 

And 

_ m/3 3/3 <R,tl^>^ +R^ <R,^> -(3R^\^\^ m^(3^k <R,tp>^ 
^" ~ ^ r^If ^ ~R^ ■ (1-A;f)2- ^ ^ 

As for the spatial derivatives we have 

Rx = ei- aif) < 'il),ei>, Rxx = (3.15) 



m <R,R^> 



^ o^ < ^, ^x >^ ^ <Rx^Rx^ kn^ < R, Rx >^ , . 



Thus 

. _ m <R,Rx>^ + < R,Ry >^ + < R,Rz >^ 

TTl < Rx, Rx > + < Ry, Ry > + < Rzi Rz > 

R^ 1-kf ^ 

km'^ <R,Rx>'^ + < R, Ry>^ + < R, Rz >^ 
~R6 (1-A;f)2 

Substitute the value of i?^ to get 



(3.18) 



v2ir/:i4 



m /?2_2c^<^^/j>2 +Q,2|^|2 ^^^^>2 ^ 3-2a|V'P + aV, 
■' ~ ^R^ 1-A;f ^ 1-A;f ^ 

A;m2 i?2_2cK^,;R>2+(;^2|^|2^^^;^>2 
:R^ (1-A:f)2 ^^^ ^^ 



Thus the second order Lh.s. of the equation is 

mP 3P<R,^>^ +R'^ <R,^> -(3R^\^W m^f3^k <R,^>^ 



R^ 






1- 


' '^R 




1 




i?6 




[i-k^y 




m 


R^ 


-2a ■ 


<Jj,R 


>^+a^i^\ 


'<^: 


R>^ 


+ 


m 
R^ 


3 


- 2a ^ 2 + a' 


2|^^ 






] 


— k^ 
- '^R 








km^ 


R'- 


-2a < 


i^,R> 


2+a2|^2 


<^,R>^ 












R^ 






(1- 


-k^y 








m(3 

i?3 


<R, 

1- 


k^ 


'^ ^n2 

«=<" 


+ 2a-a' 


'■^-). 


(3< 


R 


1- 


2 _ 
'^R 


-^^1?^)^ 




m?k 


-R' 


' + {f3' 


+ 2a- 


-a^\ip^)< 


R,i 


>2 












R^ 






(1- 


kW 

















(3.20) 



Using the relation (3 = a \il)\ — 2a we obtain 



,21 







^^ = is3 ■ 


1-fc 


m 
R 


i?4 (1 - fcf ) 


2 




As for the quadratic r.h.s- 














V'^'lah 


m^/3' 

i?6 


<R,tp>^ 


i?6 


(l-fcf)2 


+ <R,R, >2 


{1-kfy 






m2/32 

i?6 


<R,^>^ 


HI? 
R^ 


i?2- 


2a<i),R>^ 

(1- 


+a2|^ 
fcf)2 


2<^,:r>2 




{1-k^y 






9 

AM 


1 

l-^f)2 + 


rn\n2 


+ 2a- 


-a^|^^)<^ 


,^>2 

A:f)2 




Thus 




















^'^V,a/,. = - 


2 

m 


1 








i?4(l 


-fcf)2 





(3.21) 



(3.22) 



(3.23) 



Thus the field equation ( p.l| ) results in 

mP < R,ijj > 

This means i/; = 0. So the one point singularity moves with constant velocity. 
Consider now a system of n singular points. In this case the ansatz is 



(3.24) 






where 

(i)R =(r- (*V) - (*)« (»^ < ^^, (f - (*)r) > ) - «/3»^. (3.26) 

The hnear part of the equation is 



Calculate the nonlinear part 
Thus 






(/*)' = E TTTk ■ -rrfk^ < ''^' '^^ >< ''^' ''^ > • ^3-29) 



As for the spatial derivatives 






n w^ 



(')m 0)t 



<v/,v/> = E7^ffc-7^^(<^^^^-^^^^><"^^-"^^> + 

< ^"^Ry, ^'^R >< ^'%, ^^^R > + < «:R„ ^'^R >< ^-^^R,, (^'^i? > \3.3i) 
Using the relation 

^^R^ = ei - (*)a(^)^ < (*)^, ei > (3.32) 

and writing 

< ^^R^, «:R >=< ei, (')i? > -«a < «^, ^'^R >< «^, ei > (3.33) 

we obtain 

( < «:R^, ^^R >< ^-^^R^, ^^^R > ) = (^ < ei, '^^R >< ei, (■'■):R > ) - 
(■'■^a^ < ei, W;R >< (^■)^, '■^^R >< (■'■)^, ei > ) - 
(^)a( < ei, (^■):R >< (*)^, (^):R >< «^, ei > ) + 
('^a^^'^a < (*) r^, (^):R >< (')^,ei >< ^^^,^^^R>< ^^,ei > V (3.34) 



Thus the brackets in ( 3.31 ) are 

f ^ = < «:r, (j);r > -(j)a (^ < (■'■)^, »:r >< (j)^, (j^r > j - 
(')«( < (*)^, (^■):r >< (')^, «:r > ) + 

^'^a^^^a( < ^^, '^^R >< ^^, ^^^ >< (■'■)$, (^■):R > ) . 



(3.35) 



Thus, the r.h.s. of the field equation is 



< ^'^t, ^^^R >< ^'¥, ^'^R > (^"^P'^^^P + ^^a + ^^^a - '^^a^^^a < ^^, ^^^ > 



Extracting in this expression the nonhnear part we obtain 



(i)j^2^ 



(i)r 



U)r 



% l-"- ^l.i)R) i=ij '- '^ii)R ^ '^U)R 

< ^^, '^^^R >< «^, ^'^R > i^^p^^'^p + (*)« + (^'^a - '^^a^^'^a < ^^, ^^^ > ) 



Thus the field equation takes the form 

2^^^^ r.. = -k 



TiTRi 






j=l 



Wi?3 



1-A; 






(J)r 






< (')i? (^'^i? > + 



< (*V, (^'^i? >< ('V, ^'^R > (^'^P'^'^P + (')« + (^'^a - «a(^')a < ^'V, ^^'V > ) ) (3.36) 



For the approximation of slow motions it is 









0)r 



0)ij3 



4 = 1 



l-fc 



(0 



«R 



«7^i 






(j)t 



< ^^R, ^^^R > (3.37) 



The two sides of this equation are functions of an arbitrary point x. Choose the singularity 
point i = p and consider the arbitrary point x to be close to this singularity. It follows that 



^P^R -^ and ^^R -^ Rip for i ^ k, 

where Rip is a vector from the point i to the point p. 

In the l.h.s. of the equation ( p.37| ) there is one singular term 

ip)m < (P)R, (P^ > 



(3.38) 



ip)R3 1 - k 



Wr 



(3.39) 



The singular term in the r.h.s. of ( ^.37| ) is 



(i) 



m 



-k ^, y ^, < (p):r, ^^m > (3.40) 

The terms ( |3.39|) and (|3.40|) are 0{R~^ near the singularity. When these are inserted to the 
r.h.s and l.h.s. of ( |3.37D , respectively, the remainder will be 0{R~^ only if: 

< (p):r, (p)^ >= -k Y, '^'t)^ < ^''^^' ^'^^ > (3-41) 

This way, the strength of the singularity is diminished. 

Take into account that the point x is still arbitrary. Hence (|3.41|) can be valid only if 

^^^V^ = -^ETffc''^ ^3.42) 

For the limiting values in ( |3.38|) 



'"^=-*E^i^rriiw p-43) 



The second fraction differs from 1, significantly, only for small distances comparable to the 
Schwarzschild radius r = km. Thus, we neglect it. 
It follows that 



jVp ^p 



For a system of two singularity points 



-"■21 



«^ = -^-^ (3.45) 



If A; < then ( |3.44| ) and ( ^.451 ) result in attraction between the particles. The absolute value 
of k is unimportant, since it amounts to rescaling of the mass. 
This way the Newton law 



^^^mR 
j^P ^P 
is obtained. 



iP^^^J-:^^, (3.46) 
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